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ABSTRACT 

In the single-degenerate (SD) channel of a Type la supernovae (SN la) explosion, 
a main-sequence (MS) donor star survives the explosion but it is stripped of mass and 
shock heated. An essentially unavoidable consequence of mass loss during the explosion 
is that the companion must have an overextended envelope after the explosion. While 
this has been noted previously, it has not been strongly emphasized as an inevitable 
consequence. We calculate the future evolution of the companion by injecting 2 — 6 x 
10 ergs into the stellar evolution model of a 1 M donor star based on the post- 
explosion progenitors seen in simulations. We find that, due to the Kelvin- Helmholtz 
collapse of the envelope, the companion must become significantly more luminous (10- 
— 10 3 Lq) for a long period of time (10 3 — 10 4 years). The lack of such a luminous 
"leftover" star in the LMC supernova remnant SNR 0609 — 67.5 provides another piece 
of evidence against the SD scenario. We also show that none of the stars proposed as 
the survivors of the Tycho supernova, including Tycho G, could plausibly be the donor 
star. Additionally, luminous donors closer than ~ 10 Mpc should be observable with 
the Hubble Space Telescope starting ~ 2 years post-peak. Such systems include SN 
1937C, SN 1972E, SN 1986G, and SN 2011fe. Thus, the SD channel is already ruled 
out for at least two nearby SNe la and can be easily tested for a number of additional 
ones. We also discuss similar implications for the companions of core-collapse SNe. 

Subject headings: supernovae: Type la — supernovae: general 



1. Introduction 

Type la supernovae (SNe la) are among the most luminous explosive events known, and are ob- 
servable from halfway across the Universe (e.g., Rodney et al. 2012). The peak luminosities of SNe 
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la are well correlate d with the dec l ine rate of their light curves, thereby permitting their use as stan 
dardizable candles ()Phillipslll993l ; lHamuv et alJll99a : iRiess et al.lll995l ). The combinat ion of these 



two properties has made SN e la very useful for measuring cosmological parameters (jRiess et al 



19981 : iPerlmutter et al.lll999l ). As surveys have grown in scope, systematic error s are becoming the 



dominant source of unce r tainty in Type la SNe-based cosmol ogical constraints (jWood-Vasev et al 



2007 



Kessler et al 



2009 



Guv et alJl201(l : IConlev et al. 



2011 



Even as SNe la have been crucial in shaping our understanding of the universe, the nature of 
the ir progenitor system s remains theoretically ambiguous and observationally elusive (for a review 



sec 



Wang &: Hani 12012 ). While it is commonly accepted that SNe la result from the thermonu- 
clear explosion of a carbon-oxygen white dwarf (WD) in a close binary system, the nature of the 
binary companion and the chain(s) of events leading up to the SN explosion are still uncertain. 
Broadly speaking, there are two dominant progenit or models. The first is the double degenerate 



(DP) scenario where the companion is also a WD (ITutukov fc Yungelson 



1984 



1979 



Iben & Tutukov 



Webbinkl 1 1 9841 ) . The second is the single degenerate (SD) scenario where the companion is 



a non-degenerate object: a main seq uence (MS) star, a red giant (RG), a sub-giant, or a He star 



(Whelan fc Iben 



1973 



Nomotdll982 ). While the SD scenario is widely considered to be the more 



plausible channel for SNe la, there is mounting evidence suggesting that this view is incorrect. First, 
SNe la progenito r models predict outf lows during the presupernova evolution creating windblown 



cavities, however, Badenes et al 



20071 ) find that the X-ray observations of 7 young SN la remnants 



are inconsistent with this picturajj. Second, it has recently been shown that rate of WD mergers in 



suDDort of the DD scenario (Badenes &; Maoz 


2012). Additionally, giant donors seem to be ruled 


out observationally as the dominant channe 


3 (e.g.. Mattila et al. 2005: Leonard 


2007: Havden et al. 


2010: Bianco et al. 2011: Brown et al. 20121 


>; I 


A et al. 2011: Bloom et al. 2012: 


Brown et al. 2012a: 


Edwards et al. 20121: Schaefer & Pagnotta 20121"). and in this studv we do not consider them as pos- 



sible progenito r systems. An He star donor was recently ruled out for the nearby Type la SN 2011fe 
, while main sequence and sub-giant companions have mounting evidence ag ainst 



(Li et al 



2011 



them (e.g.. Beonardl 120071 : Brown et al.ll2012bl : Bloom et all 12012 : ISchaefer k Pagnottal 12012 ). but 
might still remain a significant channel for the production of SNe la. 

One a pproach to constraining the SD channel is to examine the effects of the explosion on the 



donor star ( 


Colgate 


1970: 


Cheng 


1974: 


Wheeler et al. 


19' 


75 


Frvxell k Arnett 


1981: 


Taam k Frvxelll 


1984 


Livne et al. 


1992 


Marietta et al. 


2000: Meng et al. 


20071: Pakmor et al. 


2008: 


Pan et al. 


2012Y 



While man y aspects of the interaction of a Type la explosion with the do nor star, such as its kick ve- 



locity (e.g., Ruiz-Lapuente et al. 20041), early -time shock emission (e.g., Kasen 201C : Havden et al 



2010; Bianco et al 



2011 



Brown et al.ll2012al ). the statistical properties of the surviving companion 



2 One notable exception, RCW 86, is discussed by Williams et al. 1 201lh . 

3 There are, however, a number of systems which a re considered to be candidates for the symbiotic channel (e.g. 
U Sco, RS Oph and TCrB; barthasarathv et alj|2007h . 
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population d_H an 2008; Wang &; Hanll2010l ). and the amount of material removed by the explosion 
(jMattila et al.ll2005l ; iLeonardl 120071 ). have received close attention, little has been made o f the dra- 
matic changes in the properties of the donor following the explosion. |Pan et al.l (|2012l ) recently 
emphasized that fully 65% of the mass removed from the companion is due to ablation (shock- 
mediated heat transfer) , not stripping (momentum transfer) Q so mass cannot be removed without 
heating the star because ablation and heating are not independent variables. Thus, as mass is 
ablated from the companion, the entire stellar interior i s also shock he a ted a nd the remaining en- 
velope is significantly puffed up ((Marietta et al.l |2000| ) . IPodsiadlowskil (|2003l ) modeled the future 
evolution of SN la companions by considering a constant amount of mass stripped from the donor 
star while treating the energy injected into the envelope of the comp anion as a free param eter, and 
varied it almost by an order of magnitude. While not favored by IPodsiadlowskil (|2003l ). models 
with small amounts of energy injected into the remaining star compared to the binding energy of 
the stripped material could actually become less luminous than the star's luminosity prior to the 
explosion. These low luminosity models are not, however, consistent with the simulations. 

In this paper we focus on the post-explosion state of the SD donor star. They are required 
to be more luminous than before the explosion by the physics observed in simulations of shock 
interactions. As a result, the recent observations of nearby SN remnants (SNRs) to detect or 
constrain the surviving star are far more stringent than generally assumed. In Section [2j we model 
the evolution of the donor star following the interaction with the SN ejecta and compare our results 
to previous discussions. In Section [3l we investigate nearby systems individually, first looking at 
SNRs in the Galaxy and the LMC and then discussing nearby extragalactic SNe la. Section H] 
briefly discusses the effect of shock interactions on the binary companions of core-collapse SNe. In 
Section [5] we review our conclusions. 



2. Post-Explosion Properties of the Secondary 



There have been sever a l num erical studies of the impact of the SN ejecta on a companion star 
of a SN la. iMarietta et al.l (|2000l ) performed two-dimensional simulations of the impact of the SN 
ejecta on MS stars, sub-giants, and RGs using an Eulerian hydrodynamics c ode and found tha t 
MS and sub-giant donors lost ~ 15% of their mass as a result of the impact. iMeng et al.l (|2007l ) 
then pointed out that the pre-supernova mass transfer required by the SD model would change 
the stellar structure of the companion, causing it to become more compact than an isolated MS 
star. From this they conclu ded that less mass sh ou ld be removed from the donor envelope than 
predicted by the models of Marietta et al. ( 2000 ). Pakmor et al. (|2008l ) confirmed this using a 



This result contradicts their previous work. lPan et al. (2010). which used the mixing of the SN ejecta at late-times 
in 2D simulatio ns and not trace r particles in 3D simulations to estimate the ratio of ablated-to-stripped material. 
As discussed in lPan et al.1 (|2012l ). this 2D method misclassified material that was ablated early in the simulation as 
having been stripped. 
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pact of SN ejecta 



Ivanova fc Taam 



three-dimensional smoothed particle hydrodynamics (SPH) code to model the im 
on a c ompanion whose structure was specified by the binary evolution study of 
(|2004l ). Finally, iPan et al.l (|2012l ) performed the most detailed study to date, with the inclusion of 
symmetry-breaking effects (i.e. orbital motion, rotation of the companion, and Roche-lobe overflow) 
in a simulation using the multi-dimensional adaptive mesh refinement code FLASH. They found 
that these effects again raise the mass unbound from the secondary to 16% for a MS donor. 

These studies show that there is a generic sequence of events during the impact of the SN ejecta 
on the companion. The initial impact sends a shock into the donor's stellar envelope which also 
sends a reverse shock back into the SN shell and creates a contact discontinuity between the two. 
The reverse shock develops into a bow shock around the companion which deflects much of the SN 
ejecta around the donor. Meanwhile, the forward shock penetrates the companion and the center 
of the shock slows as it encounters the steep density gradient of the stellar interior, accelerating 
once more after passing through the core. The wings of the shock, which propagate through the 
outer envelope, move faster than they do through the center because they are propagating through 
regions of lower density, causing the shock to become highly curved. The curved shock meets 
itself at the backside of the companion, depositing much of its energy in the donor's envelope and 
ablating significant amounts of mass from the surface. This creates a back pressure on the rest of 
the companion, decreasing the kick imparted by the SN blast wave. After most of the SN ejecta 
has passed the donor, the shock-heated outer envelope of the companion begins to expand and 
some of it becomes unbound. The majority of the unbound material is at low velocities and is 
mixed with the inner iron-group layers of the SN ejecta. Even though some material is unbound 
from the companion, the majority of its mass remains bound t o the star. The net he ating of the 
star is not reported, but in the standard 1 M MS model of Marietta et al.l (120001 ) the star is 
left with a hot, extended, asymmetric envelope containing ~ 10% of the mass, while the central 
pressur e is halved. There is a lmost no discussion of the subsequent appearance of the stars other 
than in Marietta et al.l (|2000j). They note that its core nuclear reaction rates will be diminished 
for a time-scale of 10 3 — 10 4 years but that the luminosity will be driven by the Kelvin-Helmholtz 
collapse of the extended envelope. Their expectation was for the star to become significantly more 
luminous (500 - 5000 L for a MS 1 M G star). 

Ideally, to study the long-term evolution of the impacted companion, one would want to model 
the three-dimensional interaction of the SN ejecta with a suite of possible donor stars until the com- 
panion is relaxed and roughly spherical. Then one would use the mean density, temperature, and 
composition profiles as the initial conditions for a one-dimensional stellar evolutio n code and follow 
the co mpanion until it settles back on the MS at the location of its current mass. I Ivanova Taam 
(|2004l ) show that t he mass of the donor s tar can range from ~ 0.6 — 1.5 Mq. We only explore the 
1 Mq MS model of iMarietta et al.l (|2000l ) because that paper contains the information needed to 
calibrate our model. Hydrodynamic simulations to calibrate our simple models at other masses or 
evolutionary states are not available, but the results presented in this paper should generally hold. 



We study this problem using the approach of 



Podsiadlowski 



(120031 ) , who followed the evolution 
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of a star which loses mass and is heated over a short perio d of time in order t o mimic the effects 
of the explosion. We used the stellar evolution code MESA (jPaxton et al.ll201ll ). allowing for fluid 



velocities. We added a module to first have a rapid phase of mass loss followed by a rapid phase 
of extra heating. In the the heating phase we also included an outer boundary R ou t, adjusting the 
mass loss rate so that any material at R > R ov x is lost in a wind on the escape time scale from 
that radius. This allowed us to treat mass loss as in Podsiadlowski (j2003l ) or to allow the heating 
phase to be responsible for some or all of the mass loss when the amount of energy added becomes 
large. Thus, the models had 4 parameters: the amount of mass lost in the first phase (AM), the 
amount of energy added in the second phase (-^heat)' the truncation radius (-R ou t)> and the time 
scale used for both the mass loss and heating phases (At). In each interval, the mass loss and 
heating rates were linearly increased and then symmetrically decreased. The internal energy added 
per unit mass was distributed as the ratio of the enclosed mass at any radius to the total mass. 
Thus more energy is dep osited at the edge of the star than near the center. This leads to results 
more consistent with the iMarietta et al.l (|2000l ) simulations than, for example, making the energy 
added per unit mass independent of radius. It also makes sense physically, as the shock will deposit 
more energy per unit mass in the lower density regions of the star. 



We first considered models where AM = 0.15 M to match the lMarietta et al.l (|2000l ) hydrogen 
cataclysmic variable (HCV) models and £ heat = 0, 0.5, 1, 2, 3, 4, 5 or 6 x 10 47 ergs. We used 
time scales of At = 1, 3 and 10 years, and -R ou t = 10, 30 and 100 R©. The £?heat = 5 and 
6 x 10 47 ergs models did not converge for At = 1 year and R ou ^ = 100 R Q due to the development 
of an oscillation as the transient ends and the outer envelope collapses. Figure Q] shows the evolution 
of the companion in luminosity and radius as a function of -E^eat anc ^ the duration At assuming 
R ou t = 100 R© after the transient phase ended. The expected luminosity and radius evolution 
depends little on At once a time period At has elapsed after the transient phase. The main 
effect appears to be an underestimate of the luminosity in that period, presumably because of the 
additional energy radiated during the transient phase. Starting with the 4 x 10 47 erg model, the 
star expands to i? ou t = 100 R©, leading to a small amount of additional mass loss. The inclusion of 
the outer boundary appears to drive the high £'j iea + models to a similar final structure so that there 
is little difference in their subsequent evolution. The luminosity of the higher -E^eat m °dels can 
be restricted by reducing R on i because it is the rapid collapse of a greatly over-expanded envelope 
that produces the highest luminosities. For example, using i? ou t = 10 R© caps the luminosity of 
the high E-^ eg ^ models at roughly 80 L©, but they sustain that luminosity for a long time and then 
roughly merge onto the evolution of the models with R ou t = 100 R©. 

Figure [2] presents the fractional change in the central temperature, density, and pressure for 
the At = 3 year, i? ou t = 100 R Q models as a function of the energy added, E-^ ea +. This figure 
shows that -^heat — 4 x 10 47 ergs is needed to reproduce the fractional change in the HCV model 



5 In the HCV model of Marietta et al. ( 2000h the binary separation was 3 R@, th e secondary was rep resented by 
a 1 Mq main sequence model, and the SN was represented by the W7 ejecta profile (jNomoto et al.llf984f ). 
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Fig. 1 . — The evolution of a 1 M MS companion after the stripping and heating phases have ended 
as a function of the energy added, ^eat (* n um ^ s of 10 47 ergs), and for a range of heating/stripping 
time scales At = 1, 3 and 10 years. These models have R ou i = 100 and the radii of the highest 
-^heat m °d e l s expand to this limit, initiating a small amount of additional mass- loss and leading 
to very similar subsequent evolution for these models. Top Panel: Evolution of the companion's 
luminosity. For the same -^heat' ^ e m °d e l s converge after a time At has elapsed. Bottom Panel: 
Evolution of the companion's radius. The ordering of the cases is the same as in the top panel. 
The peak luminosity is mostly determined by the maximum radius of the star. 
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Fig. 2. — The fractional change in the central temperature (triangles, dotted), density (open 
squares, dashed), and pressure (filled squares, solid) of the model as a function of the energy 



added E^^ . The la rge circles on each curve mark the fractional change in the HCV model of 
Marietta et alJ (2000), which is generally achieved for E-^ ea ^ ~ 4 x 10 47 ergs. Here we show the 
results for the At = 3 years, i? Q ut = 100 R© models, although the results depend little on these 
parameters. 



- 8 - 



of iMarietta et al.l (|2000l ). For the remainder of the paper we will take the E^ ea ^ = 4 x 10 ergs 



At = 3 year, and R OVL i = 100 R model as our fiducial companion model. Depositing the energy 
per unit mass unif ormly or significantly changing i^eat f rom this va hie, leads to models that do 
not reproduce the IMarietta et al.l (|2000l ) simulations well. Figure [3] shows the evolution of the 
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Fig. 3. — Similar to Figured] except showing the evolution of the companion star's effective temper- 
ature for three At = 3 years models, -^heat 1S * n un ^ s of 10 47 ergs. Note that for larger (smaller) 
amounts of energy injected into the companion's envelope, the star becomes cooler (hotter). 

effective temperature {T & q) for the fiducial model and models with higher and lower E^^- In 
general, we find that if -E^eat is large enough to initiate mass loss (the fiducial and higher E^ ea ^ 
models) then the temperature decreases, and models with insufficient heating to initiate mass loss 
(the lower -E^eat m °d e l s ) become hotter. Since the ultimate mechanism of the mass loss is ablation 
due to heating, only the models with enough energy input to initiate mass loss are realistic. The 
total energy in SNe la ejecta is ~ 1.2 x 10 51 ergs, of which ~ 3.3 x 10 ergs is incident on the star 
since, assuming the companion must overflow its Roche-lobe, a/R ~ 3. The binding energy of the 
0.15 Mq of stripped material is -E s t r ip — 6 x 10 47 ergs, so the total energy transfered to the star is 
approximately E^^ + -E^t^p — 10 48 ergs and ~ 60% of that energy is used to strip material from 
the surface, and ~ 40% heats the remaining star. Only a small fraction (/ ~ 0.03) of the incident 
shock energy is transferred to the star. Simply assuming that 3% of the shock energy is transfered to 
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the star with 60% of that energy going into stri pping mass does not perfectly reproduce the scaling 



of mass loss with a/R in Mariett a et al.l (120001 ). but it comes considerably closer to doing so than 



the analytic estimates of I Wheeler et al.l (|1975l ). More generally, the transfer of energy into heating 



the surviving star is a crucial parameter that needs to be reported for any future simulations. 

In the final analysis, the basic physics is relatively easy to understand. The shocks heat the 
star, depositing more energy per unit mass near the surface than near the center. However, the 
shocks have no physical knowledge of which material will eventually escape and the star has no 
sharp density jumps to create any features in the energy deposition. The heated material then 
expands, but with a continuous density distribution between material which ultimately escapes 
and material which remains bound. As a result, the energy deposited in material which eventually 
escapes cannot be significantly larger than the energy deposited in the material which does not 
escape. Mass is stripped beyond some effective radius -R ou t> leaving a puffed-up envelope that 
then begins to cool and collapse. The more extended this envelope, the higher the peak luminosity 
since a lower density envelope has a shorter thermal time scale. Ultimately, however, for the same 
envelope extent the star produces similar luminosities. Because there are no means in this process 
of distinguishing between e scaping and bound m aterial until the escaping material is stripped, the 



low -E^eat models following iPodsiadlowskil (|2003l ) are unphysical because they effectively make the 



heating of the unbound material independent of the heating of the bound material. When the mass 
loss is driven by ablation, the energy used to strip mass must be tightly correlated with the energy 
used to heat the surviving star and is not an independent variable. 

In this sense, alternate models where we drive the mass loss with only heating are far more 
physical. For small values of -E" nea t the star expands but does not reach R on ^ and has no mass 
loss. Larger values of -E^eat (> 4 x 10 47 ergs) expand the star to -R ut' an< ^ the stars begin to 
lose mass. The post-transient evolution of all stars that have had mass loss is very similar because 
the primary driver of their luminosity is the collapse of the envelope from R on t to its eventual 
equilibrium radius. This is illustrated in Figure HI which shows the evolution as a function of 
-^heat labeled by whether or not the stars lost any mass assuming -R ou t = 100 R Q . Ignoring the 
luminosity during the transient, much of which is radiating energy injected into the expanding 
envelope, the subsequent evolution of the mass losing transients is essentially independent of the 
deposited energy. The amount of mass lost rises steadily with E-^ eg +, and for the maximum energy 
^heat = er § s ' the mass loss is approximately the same as in the HCV model. The central 
pressure of this model is, however, too low, which means that we should concentrate more of the 
energy deposition towards the surface. 



3. Implications for Nearby Type la Supernovae 

Although challenging, searching for the donor star after a SN la explosion has significantly 
faded is a popular method for testing the SD model. This can be done either by examining sources 
near the center of historical SNRs or by observing the sites of recent SNe la in external galaxies. 
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Fig. 4. — The evolution of the luminosity for models with increasing values of -E^eat wnere the 
mass loss is drive only by the heating. The models shown by the solid lines have lost mass, while 
those shown by the dotted lines have not expanded enough to lose any mass. Note how all models 
with mass loss show similar late time evolution. The values of -E^eat are ' f rom bottom to top, 0.5, 
and 1, 2, . . . , 10 x 10 47 ergs, R out = 100 R© and At = 1 year. The vertical line marks the end of 
the heating period. The rapid collapse of the envelope leads to some initial oscillations that do not 
seem to affect the general results. 
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What has not been clearly factored into these discussions is that the surviving star must be cooler 
and more luminous than it was before the explosion. This is an unavoidable consequence of mass 
loss by ablation. Our fiducial model with £^eat — ^ x 10 47 er S s an d -^out = 100 the best match 
to simulations (shown in Figures [1] and [3|) , indicates that initially after the SN the companion has 
L ~ 300 — 1000 L and T eff ~ 5000 - 6000 K, depending slightly on At. Then, after ~ 10 years, 
we find L ~ 100 — 200 L Q and T eff ~ 5000 - 5200 K. After ~ 100 years, we see L ~ 50 - 60 L Q 
and T eff ~ 5300 - 5400 K. Finally, after ~ 1000 years L ~ 15 - 20 L and T eff ~ 5500 K. 



3.1. Supernova Remnants 



Searching for surviving stars near the centers of SNRs is a promisin g method to distinguish 



between progenitor models (jCanal et al. 



2001 



Ruiz-Lapuente et al J 12004 ). p articularly given the 



ability to determine the S N type using either X-ray spectra of the rem nant (Hughe s et al 



Badenes et al 



2006 



2008) or dust echo spectra of the actual event (e.g., iRest et al.ll2008bl ). Given 



1995 



SN la rates and SNR lifetimes, there are several examples in both the Galaxy and the Magellanic 
Clouds. The challenge is to precisely identify the expected location of the star. First, the companion 
will have a n on-negligible proper motion, which can be in any direction outwards from the SN 
explosion site ([Marietta et alj|2000l ). Second, while in general Type la SNRs are strikingly spherical, 
they will be slightly asymmetric due to differing densities of the surrounding interstellar medium 
(ISM) and/or bubbles or clumps of gas and dust in the surrounding environment, which will add 
uncertainty to the measurement of the geometric center of the SNR. Moreover, as time passes, the 
donor star fades as it returns to thermal equilibrium. 



3.1.1. SN 1 572: Tycho 's Supernova 



SN 1572, first seen o n November 6th, 1572, was well observed for two years (e.g., Tycho Brahe 



see iRuiz-Lapuente 120041 . Both the light curve (IRuiz-Lapuentd 12004 ) and the spectrum obtained 
from light echoes ( Rest et al. 2008a : Krause et al. 20081 ) show SN 1572 to have been a normal SNe 
la. In the past decad e there have been many attempts to identify a SD companion star to SN 1572. 
Ruiz-Lapuente et al.1 (|2004l ) performed the first search and found a promising candidate near the 
center of Tycho's SNR, Tycho G, a sub-giant with spectral type G2 IV, which has T q q = 5750 ±250 
K, M w 1 M , R « 1 — 3 R , and surface gravity log(j/cm s 2 )~ 3.0 — 4.0. They found that it 
is at a distance of d ~ 2.5 — 4.0 kpc, with a large radial velocity of = —108 ± 6 km s _1 , and 
a high derived proper moti on of 94 ± 27 km s _1 . However, there is significant controversy about 
the true nature of Tycho G. iFuhrmannl (|2005) propose s that Tycho G mi ght simply be a thick-disk 
star coincidentally passing in the vicinity of the SNR. Ilhara et al.l (|2007l ) classify Tycho G as a F8 
V star and argue that it could not be the donor star of SN 1572 because it does not exhibit the 
blue-shifted Fe I absorption lines expected for a star surrounded by the ejecta of a SN due to Fe I in 



the ejecta. The lack of Fe I absorption suggests that Tycho G is a foreground star. iKerzendorf et al 
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( 2009 ) measure a smaller radial ve locity (79 ± 2 km s - 1 ) and detect no proper motion, which is 
in moderate disagreement wit h the iRuiz-Lapuente et al.l (120041) > , and suggests that Tycho G is not 
associated with the SN event. iGonzalez Hernandez et al.l (|2009l ) used Keck high-resolution optical 
spectra to improve the stellar parameters, finding T g g- = 5900 ± 100 K, log(g/cms 2 ) = 3.85 ± 0.30 
dex, and ~ 80 km s _1 . They also find that Tycho G is overabundant in Ni, suggesting that it 
was polluted by the SN event. 

The rotational velocity (w ro t) °f a surviving do nor star could a lso be an important diagnos- 
tic of SNe la progenitor systems in the SD model. iMeng &; Yang] (|201ll ) predict that surviving 
companions should be rapidly rotating with v TQ ^ ~ 100 km s~ 1 because they are tidally locked 
at the time of the explosion. Tyc ho G, however, rotates slowl y with estimates of v ro ^ = 7.5 ± 2 



km s _1 and v r nt < 6.6 km s" 1 by iKerzendorf et al.l (|2009i ) and lGonzalez Hernandez et al.l (2009), 
respectively. IMeng &; Yang] (|201ll ) did not account for the interaction of the SN ejecta with the 
donor, and they suggest th i s mig ht explain the difference. In their full simulation with rotation and 
orbital motion, iPan et al.l (120121 ) found that the impact removes 48% of the angular momentum 
from a MS companion but only 18% of its mass. They estimate an equilibrium radius for the donor 
star after the SN of R ~ 2.4 Rq, which still implies a high rotational velocity of v ro ^ ~ 37 km s . 
They suggest that the donor will still be out of thermal equilibrium and thus have a larger radius. 
However, since the angular momentum scales with radius and rotational velocity as J oc M R Vj- i, 
the stellar radius must be R > 9 Rq to be consistent with the observed limits. This conservatively 
assumes solid body rotation and a self-similar radial distribution of mass in the donor. However, 
R ~ 9 R© is completely ruled out by the observed luminosity and temperature, so there remains a 
large disagreement between the observed and expected v vo i of Tycho G assuming it was the donor 
star in the SN 1572 progenitor system. 

To this list of problems we can now add that Tycho G is too low luminosity an d too high tem- 



peratu re to be the survivor of a la explosion. For any of our models consistent with lMarietta et al 
(|2000l ). the luminosity should be > 20 L Q , rather then the estimated 2 — 8 L , and the tempera- 



ture should be < 5500 K, rather than the observed 5900 ± 100 K. The only models which would 
agree with the observed values are essentially the "no heat" models that are grossly inconsistent 
with the initial conditions demanded by the simulations. Based on its luminosity, temperature and 
rotation rate, we conclude that Tycho G is not associated with SN 1572 and that SN 1572 was not 
a traditional SD SN. 



3.1.2. SNR 0509 - 67.5 



In many ways, the Large Magellanic Cloud (LMC), at a distance of 50.6±1.6 kpc (jBonanos et al 



20111 ). is a better place to look for the companions to Type la SNRs because the distances to the 
remnant and nearby stars are well-determined. Moreover, it is still close enough that Hubble Space 
Telescope (HST) or ground-based adaptive optics observations can search for lower mass main 
sequence stars without severe problems from crowding. 
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The central re gion of the LMC SNR 0509 — 67.5, created by a 1991T-like SN e la explosion 



400 ± 50 y ears ago (IHughes et al 



1995 1 : iRest et~alll2005l . [2008bl : badenes et alJboO^ . was recently 



studied by ISchaefer k, Pagnottal (|2012i ) using archival B, V, I, and Ha HST images. They were 



able to constrain the search region to l'/4 and found no ex-companion stars in the central region 
to a limiting visual magnitude of My = +8.4 mag, which they claim rules out all published SD 
models. However, this limit does not to account for t he line of sight extinctio n to the LMC 
(Ay 0.25 mag), so the ac tual limit is My = +8.15 mag ()Di Stefano fc Kilidl2012l ). Additionally, 



Schaefer &: Pagnottal (|2012l ) assumed that main sequence donor stars will be minimally affected by 



the explosion. As with SN 1572, our fiducial models calibrated to agree with iMarietta et al.1 ((2000) 
predict that the surviving star should be ~ 20 L after ~ 500 years. Crudely speaking, this means 
that a limit of My > +8.15 mag now corresponds to a limit of My > +11.25 on the unp erturbed 
progenitor star, which would correspond to a M ~ 0.25 M Q dwarf (jMarigo et al.ll2008l ). Lower 
mass companions are both rare and will be more affected by shock heating, so our conclusions 
largely hold for their scenarios as well. 



3.1.3. SNR 0519 - 69.0 



The central region of the LMC SNR 0519 - 69.0 was studied recently by Edwards et al.1 (j2012i ) 
using archival V and Ha HST images. Light echoes of SNR 0519 — 69.0 show that the SN spectrum 



is th at of a normal S Ne la that expl o ded 6 00 ± 200 years ago (jRest et al.1 120051 : lEdwards et al 



20121 ) . Unfortunately, Edwards et al.l (120121') were only a b le to constrain the location of possible 
companions to within 4'.'7. In this region, Edwards et al.l (|2012l ) found 27 MS stars brighter then 
V = 22.7 mag. The brightest star, Star 1, with V = 19.7 mag corresponds to My > +1.2, which 
would allow for a solar mass companion. However, adding the expectation that the luminosity 
must be roughly increased by a factor of 10, this really corresponds to a limit on the progenitor 
of My > +4 and a solar mass companion is only marginally allowed. The limits are much more 
strin gent for the fainter st ars in the error circle. Additionally, with a V — H a color of ~ 0, Figure 
2 of Edwards et al.l (|2012l ) shows that this star is on or close to the main sequence. Based on our 
fiducial model we would expect a ex-companion to be redward of the main sequence, suggesting 
that Star 1 is not the associated with the SNR. However, further studies of the possible companions 
with more HST observations at different wavelengths or spectra should be carried out. Currently, 
the nature of the progenitor system of SNR 0519 — 69.0 can be still be considered uncertain. 



3.2. Recent Supernovae 

Besides observing the explosion sites of historical SNRs, one could also look at the location of 
more recent SNe la once the SN has faded. This method benefits both from having very accurate 
search positions and from the greater luminosity of the surviving stars at earlier times {L ~ 50- 
— 300 Lq). The main problem is that these directly observed SN are far more distant, with the 
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closest being at a distance of > 3 Mpc. We briefly discuss four of the closest SNe la observed in 
the last century. 



3.2.1. SN 1937C, SN 1972E, and SN 1986G 



Three historic SNe la, SN 1937C, SN 1972E, and SN 1986G, are particularly good candidates 
for a study of their post-explosio n systems becaus e they are nearby and were well observed. SN 
1937C in IC 4182 at ~ 4 .4 Mpc (jSaha et all boOfih was discovered by Zwicky in the fall of 1937 
at Palom ar Observatory (Baade fc Zwickvl 119381). SN 1972E was discovered in NGC 5253 at ~ 
4.2 Mpc (jSaha et al.1 120061 ). Finally, the closest SN la discovered in the digita l imaging era, SN 



1986G, was discovere d in NGC 5128 (Centaurus A) at a distance of ~ 3.4 Mpc (jEvans et al 



1986 



Ferrarese et al.l 120071 ). Given the historical SN explosion images, it should be possible to centroid 
the SN explosion to a fraction of an arcsecond and with HST observations one could plausibly find 
or rule out a SD companion for these SNe la. 



3.2.2. SN 20 life 



The recent Type la SN 201 lf e in M 101 is the nearest la in the last 25 years. At a mere 
6.4 Mpc dshappee fc Stanekl boilh . this SN la provides an exc eptional opp o rtuni ty to put the 
strongest constraints on the progenitor system of any SNe la. iPatat et al.l (|201ll ) obtained 12 
epochs of high resolut ion spectra of SN 20 11fe, finding that it is only slightly reddened and lies in a 
"clean" environment. iHoresh et al.l ([2012]) used upper limits from both radio and X-ray obser v ation s 
of SN 201 lfe to rule out circumsteilar material from a giant companion. IChomiuk et al.l (|2012l ) 
then strengthened these radio constraints with EVLA observations, ruling out most popu l ar SD 
progenitor models. The proximity and the low extinction of SN 201 lfe allowed iLi et al.l (|201ll ) 
to place 10 — 100 times stronger constraints on the visible light from the progenitor system than 
previous studies using deep archival pre-explosion HST / ACS observations. These observations rule 
out luminous red giants (> 3.5 M©) and most helium stars as plausible donors. 



SN 201 lfe was discovered very quickly by the Palomar Transient Factory (jLaw et al.ll2009l ) on 
August 24, 2011, less than one day after explosion, allowing multi-wavelength follow-u p observations 



to be quick ly carried ou t . Wit h the early time UV, optical, and X-ray observations iNugent et al 



(|201ll ) and lBloom et al.l (120121') conclude d that the radius of the primary was smaller than 0.1 Hq 
and 0.02 R , respectively. iBloom et al.l (|2012) then went o nto sh ow th at their const r aints leave 



compact objects as the only viable primary. iBrown et al.l (|2012bl ) and IBloom et al.l (|2012l ) also 
used the very early-time UV and optical observations to constrain the shock heating which would 
occur as the SN ejecta collides with the companion. They rule out giants and main sequence 
secondaries more massive than the sun. These constraints seem to imply that the DD m odel is 
the only viable progenitor model. However, these studies both rely on the calculations of iKasen 
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(|2010l ). which assum es a Roche-l o be-oy erflowing secondary with a typical stellar structure. As 
previously discussed, iMeng et al.l (|2007l ) argued that the secondary should be smaller and more 
compact than a typical main sequence star. Such a smaller, more tightly bound companion would 
produce a weaker signal at even earlier times which could have been missed. Additionally, slightly 
less massive companions can also be too small for the SN ejecta shock to have been observed. 
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Fig. 5. — SN 2011fe light curve (points) with linear late-time decline fit (dotted black) and our 
fiducial model for a shock heated companion (dashed, see Section [2]). We note that a signature 
of the presence of a SD MS companion to a SNe la explosion would be a break in the very late- 
time light curve (~ 900 days) as the shock-heated companion becomes more luminous than the 
rapidly fading SN ejecta. We have assumed a blackbody to convert the models from $2] into Vega 
magnitudes. 



Very late time observations of SN 20 life would place further constraints on the progenitor 
system by detecting, or putting upper limits on, the existence of a shocked donor. The luminosity 
of the surviving donor star is expected to be L ~ 300 L Q , which c orresponds to V ~ 27 mag - 
In Figure E we pre sent t he BVRI light light curves of SN 20 life of iRichmond fc Smith! (|2012l ). 
Munari et al.1 (|2013l ) , and IShappee et al.l (|2012l ) along with a linear fit for the late-time decline as 
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compared to the expected evolution of our fiducial model. If the current decline rate continues, the 
shocked companion should begin to dominate the SN light in the R ~ 900 days after maximum 
light. At that point, HST observations monitoring the very late-time light curve of SN 2011fe 
should see it level off as most of the observed flux will then be coming from the shock heated 
companion and not the expanding, cooling shell of ejecta. 

However, there are other effects which might also cause the late-time SN light curve to flatten. 
First, light echoes from t he SN explo sion off nearby dust can add additional late-time flux as 



was seen from SN 1987A (ICrottsI I1988T) . Fortu nately, the environment surround SN 2011fe shows 
no evidence for gas or dust (jPatat et al.ll201ll ). Second, Auger and internal conversion electrons 
emitted in the decay of 57 Co are expected to become a significant source of hea ting ~ 750 days 
after the explosion and dominate the positron decay of 56 Co after ~ 1000 days ( Seitenzahl et al 



2009; 


Seitenzah 


1 2011 


of 


Ropke et al. 


2012 





The late-time bolometric light curve is then expected to flatten (Figure 4 



much longer than the half-life of Co (271.79 days). Thus, a single degenerate companion would 
still be observable at even later times. 



4. Core-Collapse Supernovae 

Interestingly, a similar process may also be possible during stripped core-collapse supernova 
(SN lib, lb and Ic) explosions. One of the theoretical models for producing stripped core-collapse 
super novae progenitors is using binary mass transfer to re move the outer enve lope of the primary 



(e.g., lYoon et al.1 l2010l ). a common envelope phase (e.g., iNomoto et al.lll995l ). or both. In such 
scenarios, the companion is again forced to lie relatively close to the primary, with similar issues 
about mass stripping and shock heating to the SD SNe la case. 

The maximum shock energy is essentially the same as for the la case because it is determined 
by the geometry of Roche-lobe overflow. This means that the effects on main sequence binary 
companions, are likely greatly reduced. First, the escape velocity of a higher mass star is somewhat 
larger, so less mass is stripped for a given amount of energy. Second, the internal energy of the 
higher mass star is much larger, so the fractional change in the energy created by the shock is much 
smaller. Third, the luminosity of the star is much larger, so any adjustment leading to a roughly 
constant source of additional luminosity is a much smaller fractional change in the luminosity. As a 
result, we would expect that the effects on an 8 M main sequence companion with the same shock 
energy and geometry will not very dramatic, and this is confirmed in our numerical experiments. 
Additionally, the typical ccSN binary companion (without mass transfer) can be at much larger 
distances than in a SNe la, so the effects will likely be unobservable for MS companions. 

However, if the companion has evolved into a giant star, the effects can be much greater due 
to the significantly reduced binding energy of the companion's envelope. It may even be possible, if 
the separation between the SN and a giant companion is small, for the SN ejecta to completely strip 
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the envelope of the companion, leaving only the dense core, as seen in hydro dynamic simulations 
of SNe la explosions with a giant companion (e.g., iMarietta et al.ll2000l ). The pro blem is that most 
of the time the companion star to a ccSN will still be on the main sequence (see iKochanekl 12009) 
unless the binary members were nearly exact twins (jPinsonneault Stan ek 2006). 



Conclusion 



Main-sequence (MS) donor stars surviving the explosion of a SD SNe la must be strongly 
heated and puffed up in a manner tightly coupled to their mass loss. In stellar evolution models 
th at roughly match the final state of a 1 M© star losing 0.15 Mq in the hydrodynamic simulations 
of IMarietta et al.1 ( 2000i ). the star must be heated by i% ea t — 4 x 10 47 ergs, and it commences its 
post explosion evolution with the radius and luminosity dramatically increased and the effective 
temperature modestly decreased. When mass is lost by ablation, it is not physically possible to 
avoid these effects. This means that the "leftover" stars of Tycho's SN and in the LMC supernova 
remnants SNR 0609 — 67.5 and SNR 0519 — 69.0 should still be ~ 20 times more luminous than 



before the explosion. This greatly strengthens the arguments of lSchaefer &; Pagnottal (J2012J) against 
a SD scenario for the LMC SNR 0609 — 67.5 and essentially rules out Tycho Star G as the survivor 
of SN 1572. Finally, such luminous donors should be observable with HST within ~ 10 Mpc. We 
point out that there are a number of recent nearby Type la SNe, including SN 1937C, SN 1972E, 
SN 1986G, and SN 2011fe ~ 2 years after the explosion that should be studied in detail. For the 
latter, difference imaging between late time and pre-explosion HST images should be a particularly 
powerful probe. Thus, the SD channel has been ruled out for at least two nearby SNe la far more 
strongly than assumed and can be easily tested for several more. 

Recent works have suggested the "spin-down model" as a way to avoid the existing limits by 
using rotational support of the WD to delay the explo s ion long enough to render the companion 
far fa inter than at the time of accretion (|Justhamll201ll ; iDi Stefano et al.ll201ll : iDi Stefano fc Kilic 
2012 ). In addition to several problems of fine tuning, there are no known physical systems where 
accretion leads to rotational velocities sufficiently close to breakup for rotation to be of much 
relevance to the hydr ostatic support of the sy stem. The fastest spinning WD known (RX J0648.0- 
4418, -P S pi n — 13.2 s, iMereghetti et al.1 120 111) is still spinning a t less than half the rate (-P S pi n < 6 
s) needed for rotation to be imp o rtant (jYoon &; Langerl 120051 ) . Moreover, most of the scenarios 
considered by IDi Stefano Kilid (|2012l ) lead to lower mass companions which will be even more 
drastically affected by shock heating, so our conclusions largely hold for their scenarios as well. 

Existing hydrodynamic simulations of the effects of explosions on SNe la companions have 
largely ignored the co nsequences of the explosions for the future evolution of the stars, although 
Marietta et al.1 (|2000l ) comment on the problem and provide enough information to calibrate our 
model. At a minimum, simulations should report the final outer radius of the bound material and 
the final energies of the companion star (as compared to the initial energies). Better still would be 
to report asymmetrically averaged density and energy profiles that could then be used as the basis 
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for more self-consistent models of their subsequent evolution. 
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